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A model for incomplete reconnection in sawtooth crashes is presented. The reconnection inflow 
during the crash phase of sawteeth self-consistently convects the high pressure core toward the 
reconnection site, raising the pressure gradient there. Reconnection shuts off if the diamagnetic 
drift speed at the reconnection site exceeds a threshold, which may explain incomplete reconnection. 
The relaxation of magnetic shear after reconnection stops may explain the destabilization of ideal 
interchange instabilities reported previously. Proof-of-principle two-fluid simulations confirm this 
basic picture. Predictions of the model compare favorably to data from the Mega Ampere Spherical 
Tokamak. Applications to transport modeling of sawteeth are discussed. The results should apply 
across tokamaks, including ITER. 



Sawtooth crashes in tokamaks occur when the core 
temperature rapidly drops following a slow rise [l| . Large 
sawteeth are deleterious for fusion because they spoil con- 
finement, while small sawteeth may be beneficial by lim- 
iting impurity accumulation Q. Kadomtsev suggested 
the cause is the m = l,n = 1 tearing mode [3[, where m 
and n are poloidal and toroidal mode numbers. The pre- 
dicted crash duration is the time it takes Sweet-Parker 
reconnection to process all available magnetic flux. This 
agreed with early experiments and simulations. 

Soon after, cracks in the model appeared. Crash times 
in larger and hotter tokamaks were much faster than 
Kadomtsev's prediction d, Q. Also, Kadomtsev's model 
assumes all available magnetic flux reconnects (reconnec- 
tion is "complete" ) , however experiments reveal that re- 
connection is usually incomplete [6j. Equivalently, the 
safety factor q = rB^/R^Bg does not exceed 1 every- 
where after a crash, where Rq and r are the major and 
minor radii and B v and Bg are toroidal and poloidal 
magnetic fields. 

Many models of incomplete reconnection exist, but 
there is no consensus on which, if any, is correct. Ex- 
amples include stochastic magnetic fields Q, diamag- 
netic and pressure effects at the magnetic island 0- 
Hlj ]. trapped high energy particles [pj - flij ]. a flattened 
q-profile [HI , and the presence of shear flow [H, [TtJ ■ 

The uncertainty of the cause of incomplete recon- 
nection impacts tokamak transport modeling. Low- 
dimensional transport models capture the sawtooth pe- 
riod and amplitude [l8j], but the fraction of flux recon- 
nected is an input parameter rather than self-consistently 
calculated. A self-consistent theory of incomplete recon- 
nection would improve tokamak transport models. 

In this letter, we propose a model for incomplete recon- 
nection in sawteeth due to the self-consistent dynamics 
of magnetic reconnection, building on established prop- 
erties of diamagnetic effects After describing the 
model, we present numerical simulations confirming its 
key aspects. Then, we show that the model is consis- 
tent with data from the Mega Ampere Spherical Toka- 
mak (MAST) [2(J. Finally, applications and limitations 
of the result are discussed. 

To understand why reconnection in Kadomtsev's 
model is complete, consider the m — l,n — 1 reconnec- 
tion plane sketched in Fig. [TJ The reversed (auxiliary) 
magnetic field B* is in red, the high pressure core is in 
grey, and the reconnection site is the black X. When re- 



connection begins, outflow jets (in blue) are driven by 
tension in newly reconnected field lines. Mass continuity 
induces plasma inflow from upstream (also in blue). This 
flow convects more magnetic flux (if available) towards 
the reconnection site, which reconnects. Thus, reconnec- 
tion is self-sustaining. 

We argue that the key to explaining incomplete recon- 
nection is the effect of reconnection dynamics on the pres- 
sure gradient at the reconnection site. Suppose the core 
is initially centered at the yellow X. The pressure gradi- 
ent at the reconnection site (the green arrow) is radially 
inward and relatively weak. As the reconnection inflow 
self-consistently convects the core outward, the pressure 
gradient at the reconnection site increases. The outward 
motion of the core has long been seen in observations [f| . 

In the presence of a strong out-of-plane (guide) mag- 
netic field Bh, in-plane pressure gradients lead to in- 
plane diamagnetic drifts, sketched in Fig.[T] Diamagnetic 
(w*) effects are known to stabilize linear and nonlinear 
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24 1 . It was shown [19[ that reconnection does 



2 I out ^ Vout ) 



(1) 



where v out is the reconnection outflow speed, v* Q = 
— Vp Q x YS/(q a n a B 2 ) is the diamagnetic drift velocity 
measured at the reconnection site for species a = i, e, 
and the "out" subscript refers to the outflow direction. 

We propose that the increase in v*i and v* e as the pres- 
sure gradient self-consistently increases due to reconnec- 
tion causes the left-hand side of Eq. (1) to increase. If 
Eq. (1) is never satisfied, reconnection is complete. If 
the pressure gradient becomes large enough, reconnec- 
tion ceases. Since Eq. (1) can be satisfied even when free 
magnetic energy remains, this provides a possible mech- 
anism for incomplete reconnection. This model departs 
from previous ones as it concerns pressure gradi- 

ents at the reconnection site rather than the magnetic 
islands. 

This model complements, and may explain key global 
features of, recent observations at MAST [2(j. They 
observe that |VT e | increases during a sawtooth period, 
peaking at the end of the crash (their Fig. 3), qualita- 
tively consistent with the model. They also show that 
secondary ideal-MHD instabilities are destabilized at the 
end of the crash cycle. Reconnection would also play 
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FIG. 1: (Color) Sketch of the m = 1, n = 1 reconnection 
plane. Reconnecting (auxiliary) magnetic fields B* are in red 
with the rational surface r s indicated by the dotted red line. 
Plasma inflows Vi„ and outflows v out are in blue with the re- 
connection site at the black X. The grey core moves from its 
initial position centered at the yellow X. The pressure gra- 
dient is the green arrow. The helical guide field Bh and the 
diamagnetic drift velocities v*i and v* e are shown. 



an important role in this process. When reconnection 
ceases, the electron-scale current sheet broadens, reduc- 
ing the magnetic shear in a region where |Vp| is large. 
Decreased shear is known to destabilize interchange in- 
stabilities {e.g. [HI). 

To test the model, proof-of-principle numerical sim- 
ulations are performed using F3D |2(j, a two-fluid code 
employing a two-dimensional slab geometry with periodic 
boundary conditions. This geometry is appropriate be- 
cause motion in the plane normal to the guide magnetic 
field is well described in two dimensions, toroidal effects 
are not expected to play a role on the short time scales 
in question (tens of /is), and three-dimensional toroidal 
simulations employ unphysical forcing terms to obtain 
sawteeth [27| • These simulations do not contain toroidal 
effects which lead to secondary ideal-MHD instabilities 
[2(ij because this facet of the evolution is outside the 
scope of this study. Electron pressure is evolved assum- 
ing an adiabatic ideal gas with a ratio of electron specific 
heats 7 e = 5/3. Since the relative diamagnetic speed is 
the key parameter, ions are assumed cold for simplicity. 
Magnetic fields and mass densities are normalized to ar- 
bitrary values B and po, velocities to the Alfven speed 
cao = -Bo/(47Ty0o) 1//2 , lengths to the ion inertial length 
dio = c/ujpi — (mfc 2 /AirpoZ 2 ^ 2 ) 1 / 2 , times to the ion 
cyclotron time f2^ = (Z c ffeBo/rriic)~ 1 , electric fields to 
Eq = caqBq/c, and pressures to po = B 2 /4ir, where m, 
is the ion mass, c is the speed of light, e is the proton 
charge, and Z c r is the effective atomic number. 

The coordinate system has x parallel to the inflow (ra- 
dial), y parallel to the outflow (poloidal), and z in the 
out-of-plane (toroidal) direction, invariant in the present 
two-dimensional simulations. The equilibrium has an in- 



plane magnetic field profile of a double Harris sheet, 

B J X ) = tanh ( X ~ L *I A \ _ tanh ( X + + 1, 

\ w Q J \ w J 

where L x x L y — 102.4 x 204.8 is the system size and 
Wq = 0.5 is the initial thickness of the current sheet. For 
this equilibrium, the toroidal mode number n = mani- 
festly, so the rational surfaces are x s — ±L x /4 — ±25.6. 
We focus on a single mode because there is typically a 
dominant mode in sawteeth; the n = mode is chosen for 
simplicity, but is not expected to alter the conclusions. 
The mass density is initially p = 1. The initial electron 
pressure profile is 

Pe{x) = - (pi + p 2 ) + ]z (pi - P2) x 

"tanh f* + 3 W8\_ taQh /*-3W8W 
\ w p J \ w p J 

The pressure gradient is localized near x — ±3L x /8 = 
±38.4 rather than at the rational surfaces x s . Thus, p e at 
the reconnection site is initially uniform. The length scale 
of the pressure gradient is w p — 2. The guide magnetic 
field B z (x) has a mean value of 5 with a profile that 
ensures initial pressure balance, p + B 2 /2 = constant. 

The data we present are from simulations with a grid 
scale of A = 0.05. A test simulation with A = 0.025 
confirms the resolution is sufficient. The equations em- 
ploy fourth-order diffusion with coefficient D4 = 2 x I0~ 5 
to damp noise at the grid scale; D4 has been varied to 
ensure the key physics is not sensitive to it. The electron 
to ion mass ratio is 1/25. Simulations include no resis- 
tivity because experimental crash times are faster than 
collisional reconnection times. The presented simulations 
do not employ a parallel thermal conductivity, but test 
simulations with x\\ — 0.08 reveal no significant changes. 
Tearing is initiated by a small coherent perturbation to 
the in-plane magnetic field of amplitude 0.01. It is known 
that secondary islands can spontaneously arise in recon- 
nection simulations; due to symmetry, such islands would 
stay at the original X-line [28| . To prevent this, initial 
random magnetic perturbations of magnitude 2.0 x 10~ 5 
break symmetry so secondary islands are ejected. 

The principal simulation employs p\ — 5,P2 — 25 so 
v* e will exceed v ou t when the high pressure plasma con- 
verts in. Other simulation parameters are carefully cho- 
sen: B z ^> By as is relevant to sawteeth and p e is large 
enough so the ion Larmor radius p s — c s /tt C i exceeds the 
electron skin depth d e = c/ui pe , allowing fast reconnec- 
tion to proceed [2!j[3(|. Here, c s = {^ e Z c ffT e /mi) 1 / 2 is 
the ion acoustic speed, and T e is the electron tempera- 
ture. 

Upon evolving the system, Hall reconnection occurs 
initially and the high pressure plasma convects towards 
the reconnection site as expected. The reconnection rate 
E, measured as the time rate of change of magnetic flux 
between the X-line and O-line, is plotted as the solid (red) 
line in Fig. [2ja). It increases from zero to its expected 
value near 0.1 [31( by t ~ 90, where it reaches a steady- 
state with a single X-line. (The variation between t — 40 
and 90 is due to transient secondary island formation 
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FIG. 2: (Color online) (a) Reconnection rate £ as a function 
of time t with and without a pressure gradient, (b) Diamag- 
netic drift speed v* e at the reconnection site and outflow speed 

V out VS. t. 



and coalescence.) At t ~ 195, E begins decreasing. It 
decreases to below zero, where it fluctuates for a number 
of Alfven crossing times. Thus, reconnection has shut off. 

To determine the cause, the electron diamagnetic speed 
v* e at the reconnection site is plotted as a function of 
time in Fig. ^h) as the dashed (black) line. For compar- 
ison, the outflow speed v out is plotted as the solid (red) 
line. Asymmetric outflows occur when there is a pressure 
gradient in the outflow direction [32[ , and since such gra- 
dients self-consistently generate here, v out is calculated 
as the average of the maximum electron outflow speeds 
from either side of the reconnection site, averaged over 
5d e when turbulent. 

Figure 2(b) reveals that v* e is small initially, but in- 
creases in time once the pressure gradient reaches the 
reconnection site at t ~ 140. It increases until it be- 
comes comparable to v ou t at t ~ 195 (the vertical dashed 
line), the same time E begins to decrease. Therefore, 
reconnection is throttled when Eq. (1) is first satisfied. 

To ensure diamagnetic effects occur, the out-of-plane 
current density J z near the X-line is plotted in Fig. [3] 
(a) before (t = 125) and (b) after (t — 180) the pressure 
gradient arrives, with in-plane magnetic field lines super- 
imposed. The guide field is in the — z-direction and Vp e 
is in the — x-direction. The reconnection site drifts in 
the — y-direction, the direction of v* e . Note, a secondary 
instability (recently speculated to be a drift instability 
[33}) appears. The increased variability of u* e and E af- 
ter t ~ 205 are attributed to this instability. 

To ensure the observed effect is caused by the pressure 
gradient, simulations with other pressure profiles are per- 
formed. When there is no gradient with pi = P2 = 5, 
there is no decrease in E, plotted as the dashed (blue) 
line in Fig. [2ja). The same is true for p\ = P2 = 25 (not 
plotted). When p\ = 5,p2 = 7, no drop in reconnection 
rate is observed because the maximum u* e only reaches 
~ 1, but v ou t ~ 2 so Eq. (1) is never satisfied. In sum- 





FIG. 3: (Color) Out-of-plane current density J z zoomed in 
near the X-line with magnetic field lines superimposed (a) 
before (t — 125), (b) after (t = 180), and (c) significantly 
after (t — 210) the pressure gradient reaches the reconnection 
site. The x and y axes correspond to the radial and poloidal 
directions, respectively. 



mary, the simulations confirm the basic prediction of the 
model: reconnection ceases when large enough pressure 
gradients self-consistently convect into the reconnection 
site despite the presence of free magnetic energy. 

Post-cessation features are important for the subse- 
quent dynamics. Figure 3(c) shows J z significantly af- 
ter the pressure gradient reaches the reconnection site 
(t = 210). The current layer clearly broadens as recon- 
nection stops, reducing the magnetic shear at the recon- 
nection site, as evidenced by the negative reconnection 
rate in Fig. 2(a). The reduced shear would make the sys- 
tem more prone to interchange instabilities, which were 
argued to occur in Ref. [2(| ■ 

Equation (1) provides a quantitative prediction of the 
conditions at the end of sawteeth; we assess it with data 
from MAST [20| . To transform into the plane of recon- 
nection perpendicular to the m = l,n = 1 helical direc- 
tion, the reconnecting (auxiliary) field B„ is related to 
the toroidal B v and poloidal Bg fields by 



(3I while Bu> 



(2) 



0.4 T and 



B*(r)=B e 

At MAST, Ro = 0.85 m _ 
Bg ~ 0.15 T [35j. The rational surface r s is where -B* = 
in Eq. (2), which gives r s ~ 0.32 m. This result agrees 
well with Fig. 1(a) of Ref. The helical guide field at 
r s is B h = B v (l + r s /R ) ~ 0.55 T. 

To test the model, Eq. (1) must be evaluated at the 
end of the sawtooth crash. The outflow speed scales with 
CAe, the electron Alfven speed based on the field 
upstream of the electron current layer. Assuming the 
large guide field limit with Bh S> -B* in the vicinity of 
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r s , the thickness of the electron current layer scales as 
the electron Larmor radius p e — v t h. e /^ce [11]) where 
Vth.e = {leTe/niz) 1 / 2 is the electron thermal speed and 
O ce = eB/m e c is the electron cyclotron frequency. Using 
T e ~ 500 eV at r s [H and j e = 5/3, we find p e ~ 
0.013 cm. To find £?* e , we evaluate Eq. (2) at r s ± 2p e 
.37], which gives £>» e ~5.9xlO~ 5 T, justifying the strong 
guide field assumption. Using this value gives v out 
14.2 km/s, where n e ~ 6 x 10 19 m~ 3 is estimated from 
Fig. 2 in Ref. [H. 

To estimate v* e , note | Vp e | /n e = | VT e |+T e (|Vn e |/n e )- 
The right-hand side is estimated at the end of the crash 
from Figs. 1(e), 2 and 3 of Ref. [H to be |Vp e |/n e ~ 
7400 eV/m. Then, the electron diamagnetic speed is 
V* e = \Vp e \/(qn e Bh) ~ 13.5 km/s. Equation (1) in- 
cludes ion diamagnetic effects, but complementary ion 
data is unavailable [Hj]. Assuming the ion temperature 
has a similar profile as the electrons with T e > Tj, we 
expect v* e < \v*i\ + \v* e \ < 2v* e . Thus, the two speeds 
agree rather well, showing the agreement with the data 
is also quantitative. 

As a further consistency check, we compare the speed 
of the core to the inflow speed. The the core's speed 
is estimated from Figs, l(d-f) of Ref. 20] by dividing its 
displacement (~ 0.08 m) by the elapsed time (~ 0.04 ms), 
giving a speed of ~ 2 km/s. The reconnection inflow 
speed scales like 0.1 CAi |26|, where CAi is the ion Alfven 
speed based on the field upstream of the ion current 
layer. The ion layer thickness with a large guide field 
scales like the ion Larmor radius p s [2l|. Using Z e g ~ 1 
38] and nij = 2m p for a deuterium plasma [34| . we find 
p s ~ 0.77 cm. As in the calculation of -B* e , we evaluate 
Eq. (2) at r s ± 2p s , giving B„ = 6.7 x 10~ 3 T. Then, 
CAi ~ 13 km/s, so the inflow speed is ~ 1.3 km/s. Thus, 
the inflow speed is comparable to the speed of the core, 
as predicted. 

For tokamak applications, Eq. (1) may be recast in 
terms of more familiar quantities. Assuming v ou t ~ CAe 
in Eq. (1) and rewriting Eq. (2) in terms of q and ex- 



panding to lowest order in r for a small displacement 
(2p e ) from r s , B te ~ Bgq'2p e , where the prime denotes 
a radial derivative. Thus, Eq. (1) becomes 



1 

~eB~h 



V % , Pe 



> 



2p e B e 



--q , 



(3) 



where all quantities are evaluated at r s . This expression 
is reminiscent of the condition on p' and q' for suppression 
of sawteeth derived from linear tearing theory |2ll . l39j . 

In conclusion, we have described a model for incom- 
plete reconnection in sawtooth crashes, tested the ba- 
sic physics with numerical simulations, and shown it is 
consistent with data from MAST. Interestingly, recent 
simulations of sawteeth revealed complete reconnection 
in MHD, but incomplete reconnection in extended-MHD 
with electron and ion diamagnetic effects [U Eoj ]: the 
present result may be relevant. Equation (1) may be use- 
ful for low-dimensional transport modeling, which cur- 
rently use ad hoc models to achieve incomplete reconnec- 
tion [4 If . The present results are machine independent, 
so they should apply both to existing tokamaks and fu- 
ture ones such as ITER. 

In future studies, the model should be tested with 
other extended-MHD effects such as ion diamagnetic ef- 
fects and higher x\\- The restriction on toroidal mode 
number n should be relaxed. The effect of the elec- 
tron pressure profile on the dynamics and the secondary 
(drift) instability should be addressed; this may need to 
utilize particle-in-cell simulations. Including 3D toroidal 
geometry is critical for exploring secondary ideal-MHD 
instabilities. Comparisons to multiple tokamak dis- 
charges should be done to test the scaling. 
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